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We study the properties of a capacitive 13.56 MHz discharge properties with a mixture of Ar/C2H2 
taking into account the plasmochemistry and growth of heavy hydrocarbons. A hybrid model was 
developed to combine the kinetic description for electron motion and the fluid approach for negative 
and positive ions transport and plasmochemical processes. A significant change of plasma parameters 
related to injection of 5.8% portion of acetylene in argon was observed and analyzed. We found 
that the electronegativity of the mixture is about 30%. The densities of negatively and positively 
charged heavy hydrocarbons are sufficiently large to be precursors for the formation of nanoparticles 
in the discharge volume. 
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I. INTRODUCTION 

The gas discharge in hydrocarbon mixtures is widely 
used for carbon film growth. These thin films are of great 
interest for a wide range of industrial applications due to 
their extraordinary material properties P, [3, 0] . 

The advantage of a capacitively coupled radio fre- 
quency (CCRF) discharge is that it can be used for pro- 
ducing non conducting films. Noble gases like argon and 
neon are often used as main background gases for hydro- 
carbon mixtures as their presence changes morphology of 
diamond like carbon films and leads to fewer crystalline 
defects 

Plasmochemical processes taking place in the CCRF 
discharge in a hydrocarbon mixture result in the creation 
of reactive and neutral species that in its turn leads to 
film growth on a wafer and to the formation of nanopar- 
ticles in the discharge volume. Nucleation of radicals 
due to the gas phase reactions essentially diminishes the 
rate of film growth at the substrate. However it was 
found that such particles can have their own applications 
like producing effective catalysts and composite coatings 



The gas phase reactions and nanoparticle formation 
in the capacitive 13.56 MHz discharge operating in a 
mixture of argon and acetylene were intensively stud- 
ied in recent experiments [1, [^. In these experiments 
the presence of dust particles induces periodic changes 
of the discharge properties. Detailed kinetic simulations 
of the plasma with movable dust in pure argon Jiy and 
for typical experimental conditions [ll| were performed. 
The results of these simulations allowed us to explain the 
transition between capacitive and resistive modes of dis- 
charge glow in the beginning of the dust growth cycle. In 
these simulations we used the assumption that the gas 
mixture did not change and plasmochemical processes 
were not taken into account as well as nanoparticles for- 
mation. However the inclusion of chemical reactions in 



a gas mixture can seriously impact the discharge proper- 
ties. There is a set of works devoted to the investigation 
of plasmochemical processes leading to dust formation 
for different gas mixtures (see for example with SiH4 (l^ 
and with C2H2 [13]). The main goal of this study is the 
investigation of plasmochemical processes in a 75 mTorr 
Ar/C2H2 CCRF discharge leading to dust particles for- 
mation and the influence of these processes on discharge 
properties in the initial stage of nanoparticle growth. 

The main mechanism for dust particle formation is 
considered to be the growth and agglomeration of hy- 
drocarbon chains. The role of certain radicals and ions 
(precursors) is very important, they initiate the growth 
process. Formation of precursors is usually a result of 
a reaction chain initiated by inelastic electron-molecule 
collisions (ionization, dissociation, electron attachment 
or excitation of the molecule) . The rate of a certain pro- 
cess is defined by the cross section, electron energy distri- 
bution function (EEDF) and gas pressure. Therefore the 
accurate calculation of the electron energy distribution 
function is very important for modeling plasmochemical 
processes. Calculating the EEDF is a computationally 
expensive task especially at low gas pressure when the 
EEDF is not a local function of the electrical field. 

In this work we study the plasma dynamics and plas- 
mochemical processes in a capacitive 13.56 MHz dis- 
charge in a Ar/C2H2 mixture for the conditions of 
Bochum experiments [1, . The interelectrode distance 
is 7 cm, the gas pressure is 75 mTorr, the amplitude of 
applied voltage is 92 V, the gas inlet is 8 seem for argon 
and 0.5 seem for acetylene. As a base for the chemical 
processes simulation we took the reaction set from 

The structure of the paper is as follows. The descrip- 
tion of the model and a discussion of its applicability are 
given in Sec. |TT1 The numerical results are discussed in 
Sec. mil The conclusions are given Sec. IIVI 
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II. MODEL DESCRIPTION 
A. Model overview 



space (1D3V) Particle in cell Monte Carlo collision (PIC- 
MCC) method, described in [Ti |. 



In our hybrid model the electron dynamics is described 
with the Boltzmann equation and the motion of differ- 
ent types of ions and radicals is modeled with the fluid 
approach. From the energy distribution function for elec- 
trons we calculate the generation rates of ions and rad- 
icals from the background gas. The species considered 
in our model are shown in Table HI We took the set of 
species from [l3l | for pure acetylene and extend it with 
species specific for the Ar/C2H2 mixture. Balance equa- 
tions for neutrals and ions include chemical processes. 
We take into account 146 reactions (section III HI) . This 
system of equations is solved self-consistently with the 
Poisson equation for the electric field distribution. 

TABLE I: Tlie different species included in our simulation 



Neutrals 



Ar 
C2H2 

C4H2, C(iH2 

C%H2, C10H2 
C12H2 
H2 



C2m 



Ions 
Ar+, ArH+ 
C2H+ 



Radicals 



2J 



C4H+ 



Ci2Ht 



CioHq 



C2H3 
C4H3, CfiH-j, 



'6 



Cl2-ff6 

H 

C2H 5 CaH , CqH C2H ^ C4H ^ CqH 
CsH , Cio-ff , C12H Cg-ff, CioH, CriH 



B. Electron transport 

At the gas pressure of 75 mTorr the electron mean free 
path is about 1 cm. This value is comparable with the 
width of the electrode sheath. In this case the electron 
energy distribution function is not determined by the lo- 
cal electric field. Therefore for accurate calculation of the 
EEDF we solve the kinetic equation for electron motion. 
The electron distribution function fe{t, x, v) is calculated 
from the Boltzmann equation 



dfe 

dt 



dx 



eEdh 
uip dv. 



(1) 



where Wg, mg are the velocity and mass of electrons, E is 
the electrical field, Jg is the collisional integral for elec- 
trons, which includes elastic and inelastic electron scat- 
tering with argon atoms, acetylene and other hydrocar- 
bons molecules. 

From the EEDF we can calculate the electron density 
by integration 



ne{x,t) = / feCpV 



(2) 



To solve the Boltzmann equation we use a one- 
dimensional in space and three-dimensional in velocity 



C. Electron collisions 

In an argon discharge plasma containing hydrocarbons 
the electron-neutral collisions are the main source of re- 
active radicals and ions. Since the plasmochemical reac- 
tions lead to the formation of many types of species, there 
are also a lot of electron-neutral collision types. Never- 
theless most of species have a small concentration in the 
background gas and we neglect some processes involving 
them. In fact, for modeling a mixture (5.8% of acetylene 
in the inlet) only electron collisions with Ar and C2H2 
are important for the formation of the EEDF, but we 
also take into account some collisions with molecular hy- 
drogen and with heavy hydrocarbon neutrals since these 
processes affect the chemical balance. 

We also neglect electron-neutral scattering with energy 
threshold much larger than the argon ionization energy 
since the argon ionization is the most probable inelastic 
process due to the large cross section and abundance of 
argon in the mixture. For example we neglect dissociative 
ionization of C2H2. Table HII shows the inelastic electron 
collisions taken into account in our model. 



TABLE II: the different inelastic electron-neutral collisions 
included in our simulation 





Reaction 




Type 


Rof. 


1 


Ar + e—f Ar^ + 2e 




Ionization 


[15j 


2 


C2ni?2 + e -> C2nH+ 2e, n = 1. 


.4 


Ionization 


[16j, [17] 


3 


H2+e-. H+ + 2e 




Ionization 


[18] 


4 


C2„-H'2 -f e -> C2r,H + + e, n = 


1..4 


Dissociation 


fl6] 


5 


H2+e^H + H + e 




Dissociation 


m 


6 


C2H2 -f e ^ C2H~ + H 




Dissociative 


m 








attachment 




7 




Dissociative 


m 




n — 1..6, ni + 712 = n, mi + m2 ~ 


- m 


recombination 


8 


C2/f^"="> + e -> 6-27?^"='''''' + e 




Vibrational 


[22] 


9 






excitation 


[23] 


10 -1- e ^ Ar* + e 




Electron 


[24] 


11 C2H2 + C2m + e 




level 


[25j 








excitation 





The cross sections for ionization of heavy hydrocar- 
bons taken from [l^, are calculated using the Binary- 
Encounter-Bethe model, since there is no experimental 
data concerning their cross sections in the literature. Dis- 
sociation processes do not affect the charge distribution 
directly, they only affect the chemical balance between 
hydrocarbons which participate in the ion production. 
In our simulation the dissociation cross sections of heavy 
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hydrocarbons were taken to be the same as for acetylene 
due to the lack of experimental data. 

We are interested in studying the electronegativity 
of the mixture since a large portion of Ar could make 
this mixture fully electropositive in contrast to the pure 
acetylene case. Note that the attachment process is very 
important as it is the only source of negative ions in the 
system. We also took into account dissociative recombi- 
nation of hydrocarbon cations with electrons. 

The electron elastic collisions with argon atoms [2^ . 
acetylene and heavy hydrocarbons molecules [25| were 
also considered. The hydrocarbon cross sections were 
taken to be equal to the one for acetylene). We do not 
take into account electron-electron collisions because the 
plasma density is smaller than 10^° cm^'^. The secondary 
electron emission from the electrode due to the ion bom- 
bardment is not included. Indeed, preliminary calcula- 
tions of the discharge parameters in the Ar/C2H2 mix- 
ture showed that the secondary electrons do not make 
an important contribution in the plasma balance for a 
secondary electron emission coefficient 7 — 0.1. 



D. Ion transport 

For the ions the mean free path is about 0.06 cm that 
is much less than the characteristic length of the elec- 
tric field variation. Therefore for this gas pressure the 
assumption about the local dependence of the ion energy 
distribution function on the electrical field is quite ap- 
propriate. In this case we can use the fluid approach to 
calculate the ion transport. We consider many different 
types of negative and positive ions and the use of the 
fluid approach helps to accelerate the simulations and al- 
lows us to calculate the chemical processes. Nevertheless 
below we will apply the kinetic approach to validate our 
hybrid model, (section |TTF]). 

The ion transport equations for positive ions are 

Uionn„ne - prLin^ , (3) 



dt 
duiVi 



Vi- 



dx 



eE 



dt dx 
and for negative ions 



mi 



drii driiVi 
dt dx 



Vattnnrie , 



duiVi 
dt 



+ Vi 



driiVi ^ eE 
dx rrii 



(4) 



(5) 



(6) 



where rii, Vi, rrii are the ion concentration, velocity and 
mass, n„ is the concentration of neutrals corresponding 
to the i-th ion type, Ue is the electron density, E is the 
electrical field, (3 is the dissociative recombination coef- 
ficient, Vion, Vatt and Vmom are the ionization, attach- 
ment and effective momentum transfer frequencies, re- 
spectively. In Eqs. ([3]), (O we omitted the chemical 
balance terms which will be discussed in section Hi HI 



We do not consider radiative recombination because: 
a) dissociative recombination is much faster for molecular 
ions, and b) in the case of atomic argon ions the char- 
acteristic time of radiative recombination is much larger 
than for charge exchange processes. The electrodes are 
considered fully absorptive for ions. 

In Eqs. (l6|) we should define the effective momen- 
tum transfer frequency for each ion. The following model 
was used for the description of the ion-neutral collisions. 
We take into account the elastic collisions and resonance 
charge exchange for Ar"*" ions since argon is the most 
abundant background gas. For the other types of ions 
we consider only ion-neutral elastic scattering. 

The mean ion energy distribution for Ar+ calculated 
from the kinetic model fsection fll Fp demonstrates a sig- 
nificant variation so we should take into account the en- 
ergy dependence of the collision cross section. When the 
ion energy is small, the ions interact with neutrals due 
to neutral polarization. The cross section of this process 
decreases inversely proportional to the mean relative ve- 
locity between a neutral and an ion 



(7) 



Here a is the polarization of a neutral and e is the kinetic 
energy of relative motion. 

At high ion energies the ion and a neutral interact like 
hard spheres. The cross section of this process is 



Ok 



(8) 



where ai and are the collision diameters of an ion and 
a neutral, respectively. 

The resulting elastic collision cross section is the sum 
of these two cross sections. To define effective cross sec- 
tions for all ions in the mixture we need polarizabilities 
of most abundant molecules and also collision diameters 
for these molecules and for all ions. The required pa- 
rameters for these gases are given in Table IIIIl Effective 
diameters for ions were estimated based on the diameters 
for appropriate neutrals. 



TABLE III: The Lennard- Jones parameters (ct, e) and polar- 
izations OL used in our simulation 



Molecule 


a{A) 




a/ao 


Ar 


3.54 


93.3 


11.1 


C2H2 


4.033 


231.8 


23.55 


H2 


2.827 


59.7 


5.52 


CiH2 


4.5 


320 




CeH2 


5.05 


397 




CsH2 


5.9 


470 





The collision cross sections for several ions calculated 
using this model are shown in Fig [TJ With these cross 
sections we calculate the momentum transfer frequencies 



E 



(9) 
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FIG. 1: Ion-neutral collision cross sections for different ions 
in a mixture of 5.8% acetylene and argon. 

where v is the mean relative ion-neutral velocity, is 
the concentration of the k-th component of the back- 
ground gas and <Ttr,k is the transport cross section of the 
ion-neutral collision for the k-th component of the back- 
ground gas. 

E. Poisson equation 

The Poisson equation for the electrical potential dis- 
tribution 

^cj) = A^e(^,~Yln}j, E = -^, (10) 

is solved self-consistently with Eqs. H])-® by the iter- 
ation method. The boundary conditions for the Poisson 
equation are U{d,t) = UQCos{ujt) and U{0,t) = 0, where 
Uq is the applied voltage amplitude and x=0, x=d are 
the coordinates of the electrodes. 



on the local ion energy distribution function (lEDF). In 
Fig. [5] the lEDFs calculated with the full kinetic model 
are shown at different coordinates. For our conditions the 
lEDF appeared to have an almost Maxwellian distribu- 
tion in the bulk plasma and is very different in the sheath 
region. Nevertheless we can consider the Maxwellian dis- 
tribution as a good approximation to determine the ion 
kinetic coefficients. 




E, (eV) 



F. Validation of the ion transport model 

To check the accuracy of our hybrid approach for low 
gas pressure we compared the results obtained with two 
different models. The first one is a fully kinetic model 
which includes the Boltzmann equations for electrons and 
for one effective ion type. The second one is our hybrid 
model with the same type of ions. We consider only the 
ion-neutral elastic scattering. 

The parameters of the discharge were taken the same 
as in the experiment: the interelectrode distance is 7 cm, 
the gas pressure is 75 mTorr, the frequency and ampli- 
tude of the applied voltage are 13.56 MHz and 92 V, 
respectively. 

Since the elastic scattering cross section is energy de- 
pendent, the momentum transfer frequency is dependent 



FIG. 2: Ion energy distribution function in the middle of dis- 
charge gap (x=3.5 cm) (a) and in the sheath region (x=0.75 
cm) (b). 



Fig. [3] and [3] show the ion density, ion and electron en- 
ergies obtained using the hybrid model and in the PIC- 
MCC simulations. They demonstrate a good agreement 
between the two models for the ion parameters distribu- 
tion. The electron EDF is a sensitive parameter and it 
was observed that the mean electron energy decreases in 
the case of a fully kinetic simulation. But anyway this 
deviation is about 15% which is close to the accuracy of 
a kinetic model. 
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FIG. 3: Ion density distribution obtained with the hybrid 
model (solid line) and the PIC-MCC simulations (dashed 
line) . 
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FIG. 4: Electron and ion mean energy distribution obtained 
with the hybrid model (solid line) and the PIC-MCC simula- 
tions (dashed line). 



G. Neutral transport 

We consider the following balance equations for the 
neutral density distribution 

drin 
"dt 

where Z)„ is the diffusion coefhcient, Sreactj is a source 
term which represents flow, pump and mixing processes. 
Diffusion coefficients for different species in the mixture 
are obtained from pair diffusion coefficients using Blanc 
formula 



-V(D„V(n„)) = Pnine-Vionnnne + Sreact,] , (11) 



p. 



(12) 



where Ptot is the gas pressure, Pj is the pressure of the j- 
th component of the background gas and Dnj is the pair 



diffusion coefficient. These pair diffusion coefficients can 
be obtained from Lennard- Jones parameters of neutral 
molecules using Chapman-Enskog theory [2^, [l^. For 
two species with masses rrii and rrij, and Lennard- Jones 
parameters (cr^, Si) and {aj, ej), the binary diffusion co- 
efficient at the pressure P and temperature T is given 
by 



A,, = 



3 ksT 2TTkBT/mij 



where m 



16 P Tral^nDi'f) ' 
mirrij / [rrii + rrij) is the reduced mass 



(13) 



(14) 



is the binary collision diameter and — T / Y^eleJ- ^-d(^) 
is calculated from the expression [27| 



c 



E 



G 



(15) 



with A = 1.06036,5 = 0.15610, C = 0.19300, £> = 
0.47635, S = 1.03587, = 1.52996,6 = 1.76474, i7 ^ 
3.8941 1. The required parameters for different neutrals 
were estimated based on the values given in Table IIIII 



^react,j 



can be presented as a sum 



eact^j 



s 



flowj 



(16) 



where Sfiow,j represents the gas inlet, Spump,j is the 
pumping source term and Smixj is the mixing source 
term. We use the same model as in [11] and the assump- 
tion of a perfectly stirred reactor. It gives the following 
expression 



5. 



pump, J 



T 



(17) 



where r is the average residence time of all neutrals. The 
T is adjusted in such a way that the pressure in the dis- 
charge equals the desired pressure. 



{Sflc 



;) 



Vr, 



isch 



V, 



disch 



(18) 



where Vreact is the total volume of the plasma reactor 
and Vdisch is the discharge volume. 

The electrodes are considered as fully reflective for 
molecules. For radicals we adopted a sticking model de- 
scribed in with sticking coefficients taken from p^ . 
We do not take into account wall reactions. 



H. Chemical processes 

Besides the few ions and radicals which are generated 
by electron-molecule collisions all other species appear 
as a result of chemical processes. As a base for model- 
ing chemical processes we used the set of reactions pro- 
posed in [l3| where a detailed description of the model 
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and appropriate references can be found. The further 
development of the reaction set is presented in [28J • The 
considered reactions are shown in Table HVl Note that, 
the chemical balance terms are included in the transport 
equations ©, jH), pT|) . 

TABLE IV: The chemical reactions used in our simulation 
Reaction 

Cluster growth through hydrocarbon cations 

C2H+ +H2^ C2H+ + H 

C2H+ + C2H2 C4H+ + H 

C2H+ + C2H2 C4H+ + H2 

C2H+ + C2H2 C4H+ + H 

C4H2 + C2H2 — ^ CeH^ 

C4H2 + C2H2 — > CfiH^ 

CeH+ + C2H2 ^ CsH+ 

CeH^ + C2H2 CsHq 

CgHf + C2H2 CioHg 

C8H+ + C2H2 ^ C10H+ + H2 

CioHq + C2H2 Ci2Hg + H2 

Cluster growth through hydrocarbon anions 
CinH- + C2H2 C2n+2H- +H2, n = 1..5 

Cluster growth through C2H insertion 
C2H + H2 — > C2H2 + H 
C2H -\- H — > C2H2 

C2H + C2„//2 ^ C2n+2H2 +H,n=l..b 

Cluster growth through acetylene insertion 
C2nH + C2H2 C2n + 2H2 +H,n^ 2..5 

Neutralization reactions 

C2niH + C2ii2^m ~* Cani-ff + C2n2-f^m, 711,712 = 1. .6 

Cir^^H' + Ar+ ^ Can/f + Ar 
C2„,H- + ArH+ ^ C2nH + Ar + H 

C2ni H + H2 —>■ C2nH + H + H 

Charge exchange reactions 
H2 + C2H2 — > H2 + C2H2 
Ar+ + C2H2 -^Ar + C2H+ [29] 
Ar+ + H2-> Ar + H+ [29] 

Hydrogen insertion and hydrogen abstraction 

H + C2r.H2 ^ C2„//3, 71= 1.. 3 

H + C2nH3 C2nH2 +H2, n= 1..3 

H + C2nH C2nH2, 77 = 2..6 

H2 + C2nH ^ C2nH2 +H, n = 2..6 

Other reactions 
Ar+ + H2^ ArH+ + H \29} 
C2H + C2H3 — > C2H2 + C2H2 
C2H2 + C2H — > C4H3 
C4H2 + C2H C&H-i 
C4H3 -\- H — > C2H2 + C2H2 
CgHs + H —i C4H2 + C2H2 



The main difference with respect to the reaction set 
considered for pure acetylene in [l^ is the addition of 
argon involving charge exchange reactions [29j . Due to 



these reactions a large part of argon ions is converted to 
hydrocarbon and hydrogen ions. The argon presence can 
be considered as an additional source of positive ions. 
Also due to recombination reactions it is a sink for nega- 
tive ions. Thus even if argon ions are not the most abun- 
dant, the presence of argon as a background gas makes 
the mixture more electropositive. 



III. RESULTS AND ANALYSIS 

First we analyze how the presence of a small fraction 
of acetylene and the related chemical processes affect dis- 
charge properties. 

A. Acetylene influence 

Let us evaluate the influence of 5.8% acetylene in argon 
on the discharge properties for the case when acetylene 
is just injected in the discharge and heavy hydrocarbons 
are not formed yet. So we excluded all chemical reactions 
except neutralization from the chemical balance. 

Acetylene addition causes a decrease of the positive ion 
density. In Fig. [5] the ion density distributions are shown 
for pure argon and with 5.8% addition of acetylene. In 
the second case the density of positive ions is almost 5 
times lower. It can be explained by two reasons. The 
most important one is the change of the EEDF and the 
other one is the ion-ion recombination. Indeed, in the 
case of acetylene addition negative ions are formed in 
the plasma, which can recombine with the positive ions. 
In Fig. [S]the mean electron energy distribution is shown 
for pure argon and for the mixture. In the case when 
acetylene is present in the discharge the mean electron 
energy is lower (almost twice in the bulk plasma). This 
is explained by the redistribution of the discharge power 
between the collision processes. In Fig. [7] the power con- 
sumed by ionization and excitation processes are shown. 
It is seen than even a small concentration of acetylene 
consumes a large portion of the discharge power for ex- 
citation especially in the bulk region. The reason is that 
the vibrational excitation threshold of C2H2 is small and 
these inelastic collisions efficiently cool the electrons in 
the midplane, decreasing in this way the mean electron 
energy and suppressing ionization. 

The presence of acetylene also changes the electroneg- 
ativity of the mixture. We found that the density of 
negative ions reaches 50% with respect to the positive 
ion density even for such small content of acetylene in 
the background gas (Fig. [5]). The density of negative 
ions is considerable, although the electron impact ioniza- 
tion and charge exchange reactions lead to the creation 
of much larger quantity of positive ions. However, the 
negative ions are trapped in the bulk plasma due to the 
electrical potential distribution and can rarely reach the 
electrode. In such conditions the main sink for anions 
is recombination which is rather slow. As a result the 
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FIG. 6: Mean electron energy distributions for pure argon 
and for argon with 5.8% of acetylene. 

anions can reach sufRcient concentrations in the middle 
of the discharge. 



B. Cluster growth influence 

During cluster growth the light cations and anions con- 
vert into heavy ones. The concentration of heavy cations 
increases due to their lower mobility and small recom- 
bination rate with anions. Due to the latter reason the 
concentration of heavy anions also increases. In Fig. [8] 
the electron and ion density distributions arc shown for 
the case with hydrocarbon clusters with maximum 12 
carbon atoms and for the case in the absence of clusters 
(i.e., the initial stage, before cluster growth starts taking 
place). It is seen that the electron density distribution 
over the discharge gap changes just a little. Also the 
electron mean energy distribution practically does not 



change (less than 5 %) with cluster growth. It is related 
to the fact that the positive and negative ions densities 
increase simultaneously and, therefore, the spatial charge 
distribution changes insignificantly. 

Another important sequence of cluster growth and low 
gas pumping speed is the conversion of acetylene to heavy 
neutrals C2nH2 at much larger time scale. It finally 
leads to a significant neutrals concentration variation 
i.e., a decreasing concentration of acetylene and an in- 
crease of heavy hydrocarbons and argon concentration. 
This observation is in good agreement with experiments 
[1, Q where the authors observed significant acetylene 
monomer dilution. The resulting neutrals concentrations 
are shown in Fig. [9l The calculated acetylene concen- 
tration is almost ten times less than the initial 5.8% of 
the total concentration, as a result of cluster growth and 
deposition to the walls. This conversion makes electron 
collisions with such hydrocarbons like C4H2, C6H2 more 
important. Since we use approximate cross sections for 
the ionization of C4II2, C6H2 molecules, it makes our 
results less reliable if their density increases. 
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FIG. 8: Total (positive and negative) ion and electron density 
distributions without hydrocarbon clusters(initial stage) and 
with them. 



C. Species density distribution 



c 




FIG. 10; Density distributions of the various cations. 



As a result of the electron inelastic collisions and the 
plasmochemical processes we obtained the steady state 
ion density distribution shown in Figs. [TUlandfTT] The 



10'' 



10" 



E 10" 
o 



10" 



10' 



10" 



Ar'' 



c 




X (cm) 



FIG. 11: Density distributions of the various anions. 



FIG. 9: Densities of the various neutral molecules. 

heavy cations become abundant in the middle of the dis- 
charge due to the conversion and due to their lower mo- 
bility. Near the sheath-plasma boundary the electron 
impact ionization has a maximum. As a result the densi- 
ties of Ar+ and C2Hj' ions have two peaks and the Ar+ 
ions dominate near the sheath region. It should pro- 
vide an intense argon ion flux to the electrode. We have 
compared the calculated ion flux to the electrode with 
the mass spectrum measurements from Q. Fig. [TH 
shows the calculated normalized ion flux of different ions 



to the electrode presented as a function of the ion mass. 
In Fig. [T3] the measured [1, Q positive ion spectrum 
is shown. The experimental spectrum does not demon- 
strate a significant argon ion flux in contrast with the 
calculated one. The reason is probably that in the ex- 
periments only ions are counted with energy lying in the 
fixed energy window near a certain energy value. Due 
to the varying mobility the different types of ions have 
considerably different mean energies near the electrode. 
Therefore only a part of the ions arriving at the elec- 
trode are measured. In the case of Ar+ ions the mobility 
and consequently, the energy is rather small because in 
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our mixture the Ar"*" ions participate in the charge ex- 
change resonant colhsions and in the elastic scattering. 
If we took the Maxwelhan distribution for ions then the 
flux of ions with certain energy can deviate to a great 
extent from the total flux. Considering this we obtain 
the following ion fluxes to the electrode shown in Fig. 
[Til Now a much better agreement is reached with the 
experimental data. 

The negative ions reach a sufficient concentration but 
in contrast with the positive ions the lightest anion C2H~ 
is the dominant one, with a density of about 80% of all 
anions density. This density in the middle of the dis- 
charge is twice smaller than the cation density so we can 
conclude that this mixture has significant electronegativ- 
ity regardless of a small concentration of acetylene and 
other hydrocarbons. 

Heavy cations Ci2H(^ and anions Ci2H^ are present in 
the mixture with large concentration (about 10^ cm^"^) 
and can be considered as precursors for dust particle for- 
mation. 



IV. CONCLUSION 

We developed a hybrid model for simulations of the 
13.56 MHz discharge in a C2H2/Ar mixture at a gas pres- 
sure of 75 mTorr. This hybrid model combines a kinetic 
description for electron motion and the fluid approach 
for 6 negative and 16 different positive ions. 146 dif- 
ferent chemical reactions were taken into account. We 



consider the formation of heavy hydrocarbons up to 12 
carbon atoms. Both negatively and positively charged 
heavy hydrocarbons can be precursors for nanoparticles 
formation in the discharge volume, since their densities 
are sufficiently large (wlO'' cm^'^). The total density of 
negative ions reaches about one half of the positive ion 
density. Thus a small fraction of acetylene (5.8%) in 
the argon discharge makes the mixture electronegative, 
because the negative ions are trapped in the quasineu- 
tral plasma. We have also found that injection of 5.8% 
acetylene in argon decreases the plasma density by a 
factor of 5. With acetylene added a large part of the 
discharge power is transferred into electron excitation of 
C2H2 molecules. The cluster growth does not affect the 
electron density and the mean energy to a large extent, 
but the densities of positive and negative ions increase 
since the heavy ions have smaller mobility. For the con- 
ditions of the experiments Q the important sequence 
of cluster growth and low gas pumping is a significant de- 
crease of acetylene and an increase of heavy hydrocarbon 
and argon concentrations. 
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FIG. 12: Positive ion spectrum (normalized flux of different 
cations to electrode). 
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FIG. 13: Experimental positive ion spectrum. 
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FIG. 14: Positive ion spectrum (only ions with energy near 
lOeV were taken). 



